BEC "level" for measuring small forces 
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We propose a device that consists of a trapped two-component phase- separated Bose-Einstein 
condensate to measure small forces and map weak potential energy landscapes. The resolution as 
well as the measurement precision of this device can be set dynamically, allowing measurements at 
multiple scales. 
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The mesoscopic quantum coherent matter systems, 
such as dilute gas Bose-Einstein condensates (EEC's) and 
fermion superfluids, created in cold atom traps are un- 
usually sensitive and clean many-body systems that can 
be isolated nearly completely from the influence of en- 
vironment, and possess unusual control knobs such as 
the ability to vary the inter-particle interaction at will 

. This has lead to several suggestions for BEC-probes 
that would measure, for example, Casimir forces mag- 
netic fields [3], or the earth's rotation [4]. In this letter, 
we discuss the prospect of creating a mesoscopic sized 
BEC-object (a BEC "bubble") that is free or nearly free 
- the effective potential energy experienced by the bub- 
ble's center-of-mass position is nearly constant in space. 
Therefore, the bubble's center-of-mass position becomes 
an ultra-sensitive measure of any weak external force on 
the bubble. At the same time, if the boundary region of 
the bubble is 'sharp' then its edge maps out potential en- 
ergy contours with exquisite accuracy. The freely floating 
bubble also provides an intriguing template for probing 
mesoscopic quantum behavior. By crafting shallow exter- 
nal potentials superimposed upon the flat effective 'bub- 
ble' potential, experimentalists could demonstrate many- 
body quantum interference, observe quantum Brownian 
motion and study many-body quantum tunneling. 

The actual bubble we propose to use is a trapped 
phase-separated BEC of 'B' atoms floating within a si- 
multaneously trapped, larger BEC of 'A' atoms 0, Q as 
shown in the schematic of Fig. [T] The extent of the phase 
separated BEC-B bubble is confined by the BEC-A sea 
in which it floats. By carefully choosing the trapping 
potentials and interaction strengths, the trapping, buoy- 
ancy, and surface tension forces on the bubble can can- 
cel, resulting in a freely floating mesoscopic BEC object 
near the middle of the trap. An additional increase of the 
inter-species interaction strength decreases the size of the 
interface boundary region between the BEC's, sharpen- 
ing the 'edge' of the bubble. Now the bubble's position 
is extraordinarily sensitive to the magnitude of a small 
perturbing external force. Thus, the position of the bub- 
ble becomes a measure of the magnitude of this force. 
Moreover, experimental methods based on imaging the 



shadow of the bubble as depicted in Fig.[l] can accurately 
measure the bubble displacement. Such techniques have 
proved to be yield exquisitely precise position measure- 
ments elsewhere 7]. 

Conceptually, the proposed weak force detector resem- 
bles a "level" - a device that determines the horizontal di- 
rection by detecting the cancellation of the gravitational 
force component in the direction in which a bubble in a 
fluid-filled tube can move. At the point of cancellation, 
the bubble remains stationary in the middle of the tube. 
If the effective potential of the phase separated BEC bub- 
ble is exactly flat, then the actual ground state of the 
system is a superposition of states in which the bubble 
is located in different positions. If the bubble behaves as 
a macroscopic object, then once its position is measured, 
it remains localized, 'pinned' by interactions with the 
environment - an essential ingredient of the breaking of 




FIG. 1: Schematic of the proposed BEC "level" device. The 
bubble of BEC-B (represented by a sphere) is confined within 
the BEC-A sea (shown by the ellipsoid). The arrow along 
the axis of the ellipsoid is the direction of the external force 
that causes the bubble to be displaced from the center of 
the ellipsoid. This displacement can be measured by a laser 
configuration as shown in the figure. 
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translational symmetry. If the environment interactions 
are sufficiently weak and in the presence of a shallow dou- 
ble well potential, the bubble can behave as a quantum 
object, tunneling between the wells (the bubble's surface 
tension may prevent single particle tunneling) or inter- 
fering with itself if it can follow different paths to reach 
the same final position. Such experiments that reveal the 
quantum properties of the 'bubble' can be carried out on 
a 'chip' [81] (near a surface). The main source of deco- 
herence would be the excitation of phonon modes in the 
surrounding BEC-sea. 

For the discussion in this letter, we won't require the 
extreme level of control necessary for the observation of 
mesoscopic quantum behavior. Instead, we assume that 
external interactions are present (though small) and that 
the BEC 'bubble' acts as a classical object. Even in that 
case, the bubbles can probe genuine quantum effects, for 
example, the position of two bubbles each trapped in one 
well of a double well potential would be modified from 
their equilibrium position if they experience a mutual at- 
traction, induced by Casimir-like force arising from quan- 
tum fluctuations of the surrounding BEC [9| . 

We illustrate the cancellation of bubble forces in the 
case of a 'small bubble of vanishing surface tension' - the 
simplest limit. In that case, the size of the bubble is suf- 
ficiently large to neglect the surface energy, but remains 
small compared to the length scale on which the trap- 
ping potentials VA(r) and VB(r) experienced by the 'A' 
and 'B' atoms respectively, vary. Hence, we can replace 
VB(r) « Vb(R) where R denotes the bubble's center-of- 
mass position. We also assume that the BEC gas is dilute 
and that we can describe the inter-particle interactions 
by the customary contact potentials: i'A(_B)(r — r') — 
-^A(s)^(r — r') for the mutual interactions of like bosons 
'A' (or 'B') and WyiB(r — x) = A5(r — x) for the mutual in- 
teractions of unlike bosons. In a three-dimensional BEC 
gas, the contact interactions lead to local pressures within 
the single condensate regions of average density p^(B)(r) 
equal to [Xa{b) /'^]PA(B)i^) ■ If we neglect the surface ten- 
sion of the bubble, then the condition of equilibrium re- 
quires the inside and outside bubble pressures to cancel, 
i.e. P = Aa[pa(R)]V2 = Ab[pb(R)]V2, where de- 
notes the density of an all 'A' BEC of the same chemical 
potential as the actual 'A' BEC which surrounds the B- 
bubble. Hence, the inside 'bubble' density is equal to 
/9b(R) = v^Aa/A_bPa(R-)- We write the ground state 
energy of the trapped BEC 'A' which has the 'B' bub- 
ble of Nb particles immersed as the energy Ea of the 
all 'A' BEC that has the same chemical potential as the 
one with the bubble, and an integral over the 'B' bubble 
volume f2_B centered on R. Consistent with the assump- 
tion of slowly varying trapping potentials and the neglect 
of the bubble surface energy, we approximate the energy 
densities as Va(b) (r)PA(B) (r) + [AA(i3) /^]p\b) i^) ^^^^ we 



obtain 

E = Ea+ f [VB{r)pB{r) + XBpl{r)/2 
-[VA(r)p^(r) + A^75i(r)/2]dr 
« Ea + Nb [VbCR.) - ^\a/XbVb{R) 

consistent with an effective potential energy per B- 
particle equal to Vcs,b{'R-} — Vb(R) — \/Ab/AaVa(R). 
Now, if the trapping potentials for A and B particles are 
simply proportional, VA(r) = aVBi'r), and we carefully 
tune the proportionality constant a to a — > ^'Xa/Xb 
then the effective potential as seen by the bubble van- 
ishes and it floats freely inside the BEC-A. 

This value, a = \/ Xa/ Xb is also the proportion- 
ality constant at which the bubble's position becomes 
highly sensitive to any small additional external force. 
We illustrate this for the case of a harmonic trapping 
potentials, V^B,trap(r) = [1/2)[K^x'^ + Kyy'^ + K^z'^], 
VA,trap(r) = aVB,trap(r), to which a differential exter- 
nal force F = 7z (i.e. difference of force experienced by 
the A- and B- atoms) adds an approximately linear term 
to the B-potential Vb{v) = VB,tiap(r) — 7^, so that 



Kff,B(r) = (1 - a^JXB/XA] VB,trap 



7z 



(1) 



[(1 - a^XB/XA)/2] X 



[K^x" + Kyv"" -f K,{z ~ A)2 - i^,A2] (?) 

centered around a position that is displaced along the 
z-axis from the former trap middle by a distance A — 
^[Kz ^1 — a^y Xb / Xa^]~^ ■ Hence, in equilibrium, the 
measurement A provides a direct measurement of the 
magnitude of the force, |F| =7. 

The accuracy with which the bubble's position can be 
measured depends on the sharpness of its edge. In the 
edge region if we assume the respective BEC densities to 
vary as pB(e) = ps(R)(^ - e-)/S and ^^(e) = pA(R)e/i5, 
where < e < (5 is normal to the interface area A, 
then the resulting kinetic energy contribution can be ap- 
proximated by i?kin,i5 = 4:AP[l\ + Ib]/S, where Ia{b) ~ 
h/ ^y 4:171 A{B)n'A{B)XA(B) is the coherence length of BEC- 
A(B). We also find the surface interaction energy to be 
approximately equal to Eint,s = AP{X — y/ XaXb)S /3. 
Thus the equilibrium 'edge' width obtained by minimiz- 
ing this additional surface energy E^dge = EM-a,s-\-Eiui.,s is 

given by 5 = ^ I2{t\ + /|)/(A/VA^ - 1) i. There- 
fore the 'edge' sharpness can be increased by increas- 
ing the interaction strength between unlike bosons by 
means of an interspecies scattering resonance (T^]. This 
increases the importance of the surface energy, requir- 
ing either a description that includes surface tension ex- 
plicitly, or a numerical solution of the coupled Gross- 
Pitaevskii (CP) like equations. Given the lack of space, 
we choose the latter option and solve the numerical equa- 
tions by evolving the GP-dynamics in imaginary time. 
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FIG. 2: Figure comparing an analytical calculation that ne- 
glects the surface tension and the numerical GP solution. In 
(a) the half axial length, L, of the bubble is plotted as a func- 
tion of iVi for iV2 = 2 X 10'', a = 0.9, C = 0.01, /3 = 0.9, and 
ai = 0.001, (b) the density profile obtained numerically for 
A'^i — 10^ is shown, and (c) shows the magnified view of the 
interface region in (b) . The green line indicates the boundary 
between BEC-A and -B obtained in (a). 



Motivated by recent experiments with elongated EEC's 
[l^ created by focusing a singe laser beam, we consider a 
quasi-lD configuration with Uz and ioj_, the trapping fre- 
quency in the axial and transverse direction respectively, 
and ^ = lOz/lo± « 1 for both species of bosons. 

The wavef unctions, (pA and (j)B of the two interact- 
ing condensates evolve according to the coupled Gross- 
Pitaevskii (GP) equations 



i(t)A{B) = 



-V2 



- ^^A{B) + Kl(B)(r) 



+ ^A{B)N A(B)\4'A(B) P 0A(B) + ^Nb(A) \4> B (A)\^ 4> A{B) 

where the mass of 'A' and 'B' atoms is assumed to be 
same, Va{v) = ^^^(r) = aix"^ /2(^ + y"^ j^i^ -f z^jl), all 
the physical quantities are scaled in axial harmonic os- 
cillator units, and the condensate wavefunctions 
are normalized to unity. We work in the phase sepa- 
rated regime by using X/^Xa^b = l-l- The density 
profile obtained numerically is plotted in Fig. [2l In the 
same figure we also plot the value of the half axial length, 
L, obtained by neglecting surface tension and taking the 
Thomas-Fermi approximation. For treating the quasi-lD 
situation considered here, we use a modified local den- 
sity ansatz given by ^^A(B){z) = fJ-A(B) - Va(b){z) and 

f^A(B){z) = + 4:aA(B)PA{B){z) where aA{B) is 

the two-body s-wave scattering length for atoms 'A'('B') 
ll| which determines the interaction strength by the re- 



lation A 



A{B) = 47r?i aAfB)/''^- Figure [21(c) shows excel- 
lent agreement between L and the numerically obtained 



FIG. 3: Density of BEC-A is shown as a surface plot for 
7 — 0.1. The subplots top, middle, and bottom correspond 
to q^As/Aa = 0.90, 0.93 and 0.95 respectively. The BEC-B 
bubble is located within BEC-A where pA{r) — 0. 



location of the bubble edge. 

We demonstrate the working of the EEC "level" by 
considering an external force in the axial direction F — 
7z. As argued before, we tune the inter-atomic inter- 
actions such that a^y Xb / Xa 1^ and A > \/XaXb. 
The corresponding ground state density profile obtained 
numerically is shown in Fig. [3) Here we see that the 
BEC-B bubble is displaced from the center of the trap 
by a distance A « 7/(1 — a-y/As/A^)- The boson-boson 
interaction strength as well as the trapping frequencies 
can be precisely tuned in current cold-atom experimental 
setups resulting in a large displacement of the bubble and 
thus allowing for the precise measurement of |F| = 7. 

While the above discussion was related to determining 
the force (gradient of the external potential), the pro- 
posed BEC "level" can also be used for mapping accurate 
potential energy contours. In fact such a measurement 
can be achieved using a single BEC as demonstrated in 
[1^ . This experiment related the observed BEC-density 
profile to the variation of the external potential. The ex- 
perimental accuracy of the potential was a few Hz with 
a spatial resolution of few microns. One limitation of 
the potential measurement is set by the magnitude of 
the chemical potential. One cannot arbitrarily decrease 
the chemical potential or the density of the BEC will 
drop below what is necessary for imaging the atoms or, 
if the inter-particle interaction is lowered, the coherence 
length of the BEC will increase to limit the spatial res- 
olution. In contrast, the phase separated BEC config- 
uration discussed in the present letter not only ensures 
large density inside the bubble allowing good imaging, 
but an edge width that can be reduced below typical 
BEC-coherence lengths allowing increased spatial reso- 
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FIG. 4: Density profile of the quasi-lD BEC-B for different 
values of Na- The topmost figure shows the trapping po- 
tential while the remaining four plots (from top to bottom) 
correspond to iV^ = 0, 1000, 1500, 2000. 

lution. We demonstrate this principle by considering an 
axially varying potential Vtrapiz) as shown in the topmost 
subplot of Fig. m In the same figure, the remaining sub- 
plots, show the density profile of the BEC-B as a function 
of the total number of 'A' atoms, Nj^. Here we see that 
as Na is increased, the equilibrium size of the BEC-B 
bubble decreases. This continues till a point is reached 
when the BEC-B bubble is small enough and is localized 
in the valleys of the potential landscape. Moreover, as 
mentioned earlier, the edge sharpness can be improved by 
modifying the inter-species interaction strength A. Thus 
we see that the proposed BEC "level" device can map po- 
tential energy contours arising from, for example, surface 
magnetic fields with exquisite precision. 

In conclusion, we proposed a device for measuring 
small forces (gradient of potential) by using a phase sep- 
arated two component BEC configuration. The working 
of this device is analogous to that of a level, so that we 
named it a "BEC-level". This device can also map out 
potential energy contours with exquisite accuracy. This 
device could, for instance, map potential energy varia- 
tions experienced by atoms near the surface of a metal 
caused by fluctuations of the surface fields. The BEC- 
level provides several knobs that can be tuned, allowing 
for measurements to be carried out at multiple levels of 
accuracy. 

While we have only exploited the bubble as a classical 
object, its possible range of explorations includes meso- 
scopic quantum behaviour such as many-body tunnel- 
ing, quantum interference, and quantum Brownian mo- 
tion. Along with the device proposed in this letter, the 



BEC-level can be used as a template for probing quantum 
many-body behaviour. 
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